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universe) with the Hagedorn phase softening the singular behavior associated with 
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1. Introduction 

Inflation [jl[ is a beautiful solution to several difficult problems in cosmology; the 
horizon problem, the flatness problem, and the monopole problem. In field theory 
however, the standard way to obtain inflation is to add a positive cosmological con- 
stant (which has a negative pressure p = —p). This ingredient is, without doubt, the 
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least attractive feature of standard inflation and it is generally extremely difficult to 
control its adverse effects (e.g. the graceful exit and moduli problems). It is worth 
asking therefore if there are other forms of matter that can have a negative pressure 
and hence give an accelerating scale factor. In this paper we introduce a candidate 
that actually has this property - open strings on D-branes at temperatures close to 
the string scale. 

At sufficiently high temperatures and densities fundamental strings enter a cu- 
rious 'long string' Hagedorn phase 0, |3|, f|, [5], || |7|]. To date applications of this 
phase have been quite limited in string cosmology [[J because the thermodynamics is 
governed by the finite temperature partition function. A rigorous analysis therefore 
requires nothing less than solving the string system in a cosmological setting, a dif- 
ficult problem that might at best be tractable only in a few special cases. Moreover, 
in order to understand the effect of macroscopic phenomena such as winding modes 
we need the microcanonical ensemble (as we shall discuss) - an ensemble that does 
not particularly lend itself to cosmological applications. 

In this paper we show that, for certain systems, it is possible to bypass these 
technical difficulties by using a classical random walk picture to model the behaviour 
of the strings in cosmological backgrounds. The particular systems we will focus on 
are D-branes in the weak coupling limit || . In this limit one can separate the energy 
momentum tensor into two components; a localized component corresponding to the 
D-brane tension, and a diffuse component that spreads into the bulk corresponding 
to open string excitations of the brane. (At the risk of causing confusion we will 
often refer to the latter as a 'bulk' component.) We will in addition allow a bulk 
cosmological constant, although our focus in this paper will be on the cosmology when 
the combined effect of the brane and bulk cosmological constants is subdominant. 

The random walk picture allows us a first glimpse at the cosmological effects of 
a primordial Hagedorn phase of open strings on branes and we find two interesting 
types of behaviour. 

• The first we call Hagedorn inflation. We will show that the transverse 'bulk' 
components of a D-brane's energy-momentum tensor can be negative. If all of 
the transverse dimensions have winding modes, this negative 'pressure' causes 
the brane to power law inflate along its length with a scale factor that varies 
as a ~ even in the absence of a nett cosmological constant (as shown in 
eq.( |5.2|) ). If there are transverse dimensions that are large (in the sense that 
the string modes are not space-filling in these directions), then we can find 
exponential inflation (as shown in eq.( |5.3| )). 

• If there is a small but negative cosmological constant, the Universe can enter 
a stable but oscillating phase; i.e. a 'bouncing' universe. The nett effect of the 
Hagedorn phase is to soften the singular behaviour associated with the collapse. 



2 



Such singularity smoothing is a familiar aspect of strings, but the nice feature 
here is that we find it in a purely perturbative regime. 

We should at this point also emphasize a general observation that we make, 
namely that the diffuse stringy component can have a dominant effect on the cosmol- 
ogy even in the weak coupling limit. At first sight this may be somewhat surprising 
given that the intrinsic tension energy of a D-brane goes like pb r ~ l/<7s ~ l/£ where 
g s is the string coupling and k is the effective gravitational coupling. However, we 
will see that the cosmological effects of the two components are proportional to k^p\ r 
and k 2 p for the brane and diffuse 'bulk' components respectively. Then the contri- 
bution of the brane component is proportional to k A p 2 r ~ k 2 . Since j->,p>l to 
be in the Hagedorn phase, if e.g. the transverse volumes are of order unity in string 
units, then the contribution of the bulk component is k > k 2 p > k 2 . Hence the 
cosmological effect of the diffuse bulk component can be dominant when g s becomes 
small. 

We begin in sections 2 and 3 by deriving the energy-momentum tensor, the prin- 
cipal ingredient of Einstein's equations. Since the results can be understood rather 
intuitively, this part of the discussion is organized so that cosmologists (and indeed 
anybody else) can skip the bulk of sections 2 and 3 concerning string thermodynam- 
ics and proceed directly to the energy-momentum tensor which is summarized at the 
end of section 3. 

The thermodynamic discussion of section 2 gives a detailed introduction to the 
behaviour of both type I and type IIA/B open strings on D-branes as calculated 
from the microcanonical ensemble in a flat background. In particular we discuss the 
importance of macroscopic modes such as winding modes. Much of this section is 
a collation of results from ref. |5[]. We then reintroduce the classical random walk 
picture paying special attention to the meaning of quantities such as average string 
length. 

In section 3 we use the thermodynamic results to calculate the energy momentum 
tensor of the Hagedorn phase. enters into the higher dimensional Einstein's 
equations and determines the cosmology, and in particular we show that open string 
winding modes gives negative transverse components. For convenience the results 
for are summarized at the end of section 3 where we also discuss the heuristic 
interpretation of this negative pressure. 

Armed with the energy-momentum tensor, we examine the resultant cosmology. 
In sections 4 and 5, we solve the equations of motion with various ansatze, and find 
the advertised inflationary behaviour as well as bouncing solutions with singularity 
smoothing behaviour. 

We will, purely for definiteness, consider adiabatic systems in solving the evolu- 
tion equations of the universe. Under the assumption of adiabaticity, the inflationary 
growth period drives down the temperature of the system; eventually the temperature 
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drop causes the universe to leave the Hagedorn regime, and consequently inflation 
ends automatically. However adiabatic systems are probably unable to provide a 
realistic scenario with sufficient inflation. In section 6 we therefore discuss how, in 
non-adiabatic systems, inflation can be sustained. We conclude in section 7. 

2. The Hagedorn phase and random walks 

This section presents some background thermodynamics needed to get T^ v (which 
enters into Einstein's equations) in the Hagedorn regime. In section 2.1 we review 
the thermodynamic properties of D-branes in toroidal compactifications. These com- 
pactifications allow us to use the microcanonical ensemble, which is defined in terms 
of global parameters such as total energy and volume. The importance of the micro- 
canonical ensemble is that it allows a rigorous understanding of the effect of winding 
modes. This understanding enables us in sections 2.2 and 2.3 to extend the analysis 
to more general universes of any shape; in particular, we can determine those cases 
in which thermodynamics and hence cosmology can be studied. Essentially, we will 
argue that thermodynamics only makes sense in a cosmological setting for those sys- 
tems in which the canonical and microcanonical ensembles agree. Where this is not 
the case the systems are dominated by large scale fluctuations. 

For the systems in which the two ensembles agree, we will derive a expression 
for the partition function based on a heuristic random walk argument. In particular 
this expression gives a geometric understanding of the partition function that allows 
us to discuss its validity in limits of high energy density, small volume, etc. 

In section 3, we then use this partition function to find the energy momentum 
tensor, T^ u . Readers whose principal interest is cosmology may wish to read subsec- 
tion 2.1 and then proceed directly to the summary of the results for T^ v in section 
3.4. 

2.1 String thermodynamics and the Hagedorn phase 

The Hagedorn phase arises in theories containing fundamental strings because they 
have a large number of internal degrees of freedom. Indeed, because of the existence 
of many oscillator modes, the density of states grows exponentially with energy e, 
u>(e) ~ e e^ H£ , where the inverse Hagedorn temperature (3h (where (3 = 1/T) and 
the exponent b depend on the particular theory in question (for example heterotic or 
type II) [[|. For type I,IIA,IIB strings the numerical value of the inverse Hagedorn 
temperature is (3h = 2^/2rr in string units. It is easy to see that thermodynamic 
quantities, such as the entropy, are liable to diverge at the Hagedorn temperature; 
obtaining the partition function Z with the canonical ensemble and multiplying by 
the usual Boltzmann factor e _/3e , one finds an integral for the partition function (for 
large e) Z oc / dee~ b e^ H ~^ £ which diverges at (5 = /3# for b < 2. 
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If an infinite amount of energy is required to reach T#, then we say that the sys- 
tem is limiting. If not then the system is said to be non-limiting. Already the simple 
canonical ensemble above indicates that the Hagedorn temperature is a limiting tem- 
perature for the b < 2 cases [[|. The remaining systems seem to be non-limiting, 
and until recently it was thought that this might imply some kind of Hagedorn phase 
transition (drawing strongly on the analogy with the quark-hadron phase transition) 
to more fundamental degrees of freedom. However, in ref.|| it was shown that all 
string systems are in fact limiting, including arbitrary Dp-branes (i.e. open strings 
attached to p dimensional defects) . 

The limiting/non-limiting question is a rather subtle one, and in order to resolve 
it one must first allow for large scale fluctuations by using the micro canonical ensem- 
ble (as pointed out in the early papers of Carlitz and Frautschi |2|) and second, retain 
full volume dependence (as pointed out for heterotic strings by Deo et al f|] and for 
open strings in ref. 0). Once both of these factors are included, it becomes clear 
that in all cases the Hagedorn temperature is truly a limiting temperature rather 
than an indication of some sort of phase transition. Perhaps the best evidence for 
this is that the Hagedorn phase completes (by entropy matching) a phase diagram 
which includes other non-perturbative phases such as black holes. (There is a sense 
in which the entire Hagedorn phase can be thought of as a first-order phase tran- 
sition from Yang-Mills/supergravity degrees of freedom to black-branes/black-holes. 
We return to this point later.) 

Let us summarize the rigorous results from ref.0 for open and closed strings 
in a toroidally compactified space. The most direct route to the thermodynamics 
is to evaluate the one loop partition function with the Euclidean time coordinate, 
r, compactified with radius (3. However for the random walk discussion later, it is 
useful to begin with the density of states u(e) of a single string of energy e from 
which the same results are obtained; 

uli e ) = l' i "^ V 'y0^ ° pm (2.1) 
(fee)Tc+ 1 closed. 

In the above V± and V\\ are the volumes transverse and perpendicular to the D-brane. 
For open strings, 

7o = y - 1, (2.2) 

where d is the number of dimensions transverse to the brane in which there are no 
windings and V a is the volume of this space (if there are windings in all dimensions, 
then V Q = 0(1) in string units). Similarly, for closed strings, 

7c = § (2.3) 

where d c is the number of dimensions in which closed strings have no windings and, 
again, V c is the volume of this space. Note that 7 D and 7 C are e- dependent critical 



5 



exponents, because winding modes are quenched or activated depending on the string 
energy. Below we shall use 7 to stand for either 7 C or 7 Q as appropriate. 

We now collect the results obtained in ref.|| in the microcanonical ensemble 
working in an approximation to the thermodynamic limit. The two main types 
of behaviour are the single-long-string or non-limiting NL behaviour (see discussion 
below where we discuss the consequences of placing 'non-limiting' systems in thermal 
contact with limiting systems), with entropy density 



of the form 



onl[ 7 ] « (3 h p 



a = S/V\\ 

1 + 7 



log (p), 



(2.4) 
(2.5) 



and the various types of 'limiting' behaviour Lpy] 



P H p + 2Jp/V ± if 7 = -l 



(3hP+^P^ 



p H p + log (p) 



if 7 = - 
if 7 = 



1 1 

2' 2 



(2.6) 



(3 H p - e p if 7 = 1, 

up to positive constants of O (1) in string units. Here 

P = E/V\\ , 



[2.7) 



where E is the total energy of the system. Note that the microcanonical ensemble 
results are expressed in terms of global parameters such as total energy which are 
valid in a toroidal compactification. We shall find that it is only possible to generalize 
the discussion for those systems in which the canonical and microcanonical results 
agree. 

The L[-l] system is our 'standard' high energy regime. From eq.(2), we see 
that it corresponds to d a = 0, so that all dimensions have windings and it is the 
system which is always reached provided that the volumes are finite and the energy 
density is high enough. These systems are equivalent to D — 1 branes where D is the 
total number of space-time dimensions (in other words freely moving open strings). 
(Once there are many windings, we can T dualize the Dirichlet directions so that 
they become Neumann directions much smaller than the string scale - the winding 
modes become a spectrum of Kaluza-Klein modes indicating open strings which are 
energetic enough to probe all of the D — 1 Neumann dimensions - even the small 
ones.) 
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There are two other sorts of behaviour, marginal limiting ML and weak limiting 
WL with entropy densities 



%l ~ f3 c p + 



oml ~ (3 h p- P 



2 



2D-3 



P~ 



(Vb-i) 




(2.8) 



respectively, where / = V\\/V±. 

Whether a system is limiting or not is a function of the global parameters of 
the system such as total energy, E, and volumes V±, Vn. For example, imagine 
increasing the energy of an L[-l/2] system. These systems are characteristic of an 
intermediate energy phase of open string systems possessing one large transverse 
dimension without windings. As the energy is increased windings will eventually be 
excited in this direction. If the transverse dimension is of order R±, this happens at 
an energy threshold E = R\, and for higher energies the thermodynamic behaviour 
changes to that of an L[-l] system, the universal high-energy regime 4 . 

The NL behaviour is so called because the temperature is higher than T# as can 
easily be seen when we calculate the temperature of the subsystems from 



where f2« is the microcanonical density of states in the subsystem i. The non-limiting 
systems formally obey 



However it is important to realise that this simply means that the NL regimes are 
transients for a finite ten-dimensional volume. Equilibrium can never be achieved 
when they are in thermal contact with the surrounding (colder) limiting system of 
closed strings. 

We can compare the microcanonical ensemble results reviewed in the previous 
two pages to the canonical ensemble by taking the thermodynamic limit {i.e. letting 
V\\ become infinite whilst keeping the density on the brane p constant). If we consider 
strings attached to a single brane embedded in D = 10 space time dimensions, then 
there is universal agreement between the canonical and microcanonical ensembles 

4 Note that the decoupling of winding modes implied by these equations of state is universal and 
in particular independent of the zero-temperature vacuum contribution to the partition function 
which may or may not be finite as the transverse radius becomes infinite ||. This is because in 
the finite temperature piece of the partition function the winding modes are always accompanied 
by a Boltzmann-like suppression factor. A related point is that as well as the above contributions 
there is a suppressed contribution to the entropy density of order <tym ~ j 7 - 0- 2 7«— 3 corresponding 
to a Yang-Mills gas in the D — 2~f Q — 3 dimensions with windings and KK modes. In the Euclidean 
approach this contribution comes from a UV cut-off in the Schwinger integral and represents a field- 
theoretic contribution which is subdominant for temperatures close to 1/ (3h but again is universal. 



Si(Ei) =logfi i (£ i ) 
Tr 1 =pi = dSt/dEt 



(2.9) 



T L <T H < T NL . 



(2.10) 
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when there are < d± < 4 transverse dimensions. In this case heading towards the 
thermodynamic limit quenches winding modes - the transverse dimensions can be 
thought of as effectively infinite in this limit. However when d± > 4 (Dp-branes in ten 
dimensions with p < 5), modes as large as the transverse dimensions can be quenched 
if V\\ <C R?l, giving NL behaviour, or activated into a L[— 1] system for V\\ £ \fVf_ 
(when the saddle point approximation is valid). Taking V\\ — > oo with p fixed and 
R± constant leads to the onset of L[— 1] behaviour in the thermodynamic limit. 
WL behaviour occurs in an intermediate region of parameter space where windings 
are being quenched in an L[-l] system and the usual saddle point approximation is 
invalid. Finally closed strings in the thermodynamic limit have a critical dimension 
D — 3, where D is the number of spacetime dimensions. For D < 3, we get standard 
canonical behaviour, L[(Z) — l)/2]. On the other hand, for D > 3 winding modes 
generate 'marginally limiting' behaviour ML. 

Some examples of different situations are shown (somewhat impressionistically) 
in figure 1 where the dimension dependence is evident. The entropic preference of 
strings for branes with higher dimensionality has the obvious interpretation that, 
when the strings are volume filling, the larger dimension branes 'cut across' more 
strings. 

Note that in the microcanonical ensemble we always use the definition of tem- 
perature as derived from the microcanonical definition of entropy, (3 = OeS^E, V, N) 
where E is the energy of the total system. The discussion above illustrates the main 
reason for this rather convoluted set of definitions (Z((3) — > Q — > S(E) — > (3); the 
effect of winding modes cannot fully be explored using the canonical ensemble. For 
example the canonical ensemble simply doesn't know about the ML and WL sys- 
tems. The microcanonical ensemble also allows us to look at cases where parts of 
the system never come into equilibrium. This is true for instance for a D3-brane 
placed in a reservoir of closed strings. As we have seen, for large transverse volumes, 
a D3-brane formally has a temperature that is higher than T# and so any ther- 
mal energy that it may have quickly evaporates into the surrounding bath of closed 
strings. Only the microcanonical ensemble allows one to compare the temperatures 
of different subsystems. Indeed the failure to do so leads to the false conclusion that 
the system as a whole might be non-limiting. Moreover, strictly speaking the ther- 
modynamic limit {i.e. the assumption on which the canonical ensemble rests) does 
not actually exist for open strings. The phase diagram derived in ref. H] tells us that 
when trying to take the thermodynamic limit one inevitably encounters a black hole 
phase. Thus it makes more sense to define the system with global quantities such 
as total energy, and work in finite volumes taking an approximation to the naive 
thermodynamic limit if possible. 

2.2 The random walk interpretation 

Despite its many advantages the microcanonical ensemble has one distinct disadvan- 
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tage; in many cosmological backgrounds global properties such as total energy cannot 
be defined. We therefore need to be able to model these aspects of the thermody- 
namics with a more flexible approach. The random walk interpretation provides just 
such flexibility. For example, we can easily understand the NL/L cross over if we 
interpret a string of energy e as a random walk of length ~ e in string units. The 
size of such a random walk is yfe so that, roughly speaking, to have modes of size R± 
we require that the total energy E = pV\\ ~ Vu ~ R\. In this classical picture the 
NL/L cross over occurs precisely when there are modes as large as the transverse 
volume 5 . One helpful feature of the random walk interpretation is that we no longer 
have to bother about the precise topological properties of the space, merely its size 
in string units. On the other hand what we have gained from the rigorous analysis of 
the microcanonical ensemble is confidence in our understanding of which systems are 
stable, when one can use naive canonical ensemble results {i.e. as an approximation 
to the thermodynamic limit when the systems are limiting) and when the random 
walk interpretation is valid. 

An analysis of random walks was the approach taken by Lowe and Thorlacius 
for closed strings and Lee and Thorlacius for open strings attached to D-branes 0. 
These authors studied the Boltzmann equations for the average number of interacting 
strings sections of different lengths. The equilibrium solutions for closed string loops 
of length I for example show that the (single string) density of states must be of the 
form u)(e) ~ ^e 13 " 6 , where e = a£ and the temperature is related to the average total 
length L of string in the ensemble by 



This supports the classical interpretation of the string as a random walk with, 
e = cr£, where a is the string tension encoding information about, for example, the 
step-size. In ref. || it was shown that such a random walk interpretation also accounts 
for the volume dependence as follows. First consider the distribution function ui{e) 
for closed strings in D large space-time dimensions. The energy e of the string is 
proportional to the length of the random walk. The number of walks with a fixed 
starting point and a given length e grows exponentially as exp e). Since the walk 
must be closed, this overcounts by a factor of the volume of the walk, which we shall 
denote by 



5 We should caution against taking this picture too literally although, as we are about to see, it 
does better than might be expected in explaining the thermodynamics. In particular the way to 
measure the 'size' of the string is to put it in a box with a certain energy and see at what radius 
it begins to excite winding modes. This should be borne in mind in order to avoid getting into 
circular arguments. 



Ph = P- t - 



(2.11) 



V(walk) = W . 



(2.12) 
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Finally, there is a factor of Vd-i from the translational zero mode, and a factor of 
1/e because any point in the closed string can be a starting point. The end result is 

1 e@ H£ 

enclosed ~ Vd-1 (2.13) 

e W 

Now, the volume of the walk is proportional to e( £)_1 )/ 2 if it is well-contained in the 
volume (R 3> or roughly Vp-i if it is space-filling (R <C y/e). From here we get 
the standard result || f4[. We have 

^O)closcd ~ Vb_i (jD+1)/2 (2.14) 

in D effectively non-compact space-time dimensions, and 

^Wcd ~ — ^— (2.15) 

in an effectively compact space. Note that these results agree with those of eq.(l); 
for example to obtain the same result as eq. (|2.14 ), we take d c — D — 1, the number 
of space dimensions. 

We can customize this analysis for open strings attached to a brane by a slight 
modification of the combinatorics. (The more general case of intersecting p, q branes 
was discussed in ref. 0). The leading exponential degeneracy of a random walk of 
length e with a fixed starting point in say the Dp-brane is the same as for closed 
strings: exp(/3#e). Fixing also the end-point at a particular point of the Dp-brane 
requires the factor 1/W to cancel the overcounting, just as in the closed string case. 
Now, both end-points move freely in the part of the brane occupied by the walk. 
This gives a further degeneracy factor 

Wf (2.16) 

from the positions of the end-points. Finally, the overall translation of the walk in 
the excluded transverse volume gives a factor Vy/W||. The final result is: 

V\\ 

^OOopcn ~ 7TT ex P (Ph£)- (2.17) 
W j_ 

Thus we see the sensitivity of the density of states to the effective volume of the 
random walk in the transverse directions. If the walk is well-contained in these 
directions, (R± 3> \fe), then we find W± ~ e d ±/ 2 and 

u(e) opcn ~ exp(p H e) R± > y/e. (2.18) 

On the other hand, if it is space-filling in the transverse directions (R± \fe), the 
transverse volume of the walk is just W± ~ V± and we find 

^(£)open ~ jy- expipHc) R± (2.19) 
v± 

in agreement with eq.( |2.1| ). 
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2.3 Random walks in a cosmological background 



The above thermodynamics involved a Euclidean metric. How can we adapt the 
results for a cosmological background? First we assumed that the density of states 
increased as exp(/3#e). In a sense the parameter e is now no longer the energy 
but becomes merely something to be integrated over. In the case of a non-trivial 
metric the most natural interpretation of e is that it is the proper length of a string 
in the bulk and certainly we can always go to the local inertial frame in which a 
small portion of the string has the usual Euclidean energy = length equivalence. In 
addition it is clear that the average number of strings of a given proper length is well 
defined and so the arguments of Lee, Lowe and Thorlacius then give us the correct 
distribution of proper string lengths in equilibrium. Furthermore eq. (|2.11 ) then gives 



us a working definition of temperature in terms of average proper string lengths. 

To put these arguments on a firmer footing we now derive one aspect of the 
strings' behaviour in a cosmological background that will be useful in the discussions 
that follow; namely that in a slowly varying background a classical string of proper 
length t ~ e occupies a volume of proper-size \fe. 

Consider a classical string of proper length £ = e/a (where a ~ 1 in string 
units) with one end point fixed in a D — 1 dimensional space. From our observations 
above we would expect such a string to have a density of states u>(s) ~ e 13 " 6 . The 
crucial point is that we can arbitrarily divide this string into N small strings each 
of which has one end free and one end that is fixed to the end of the previous 
string. Consequently the density of states of the large single string is the same as 
that of N small strings each of which has one end fixed and with energies obeying 
J2i e % — £ where £j is the energy of each string in its local inertial frame. By choosing 
a sufficiently large N we may always use the flat space approximation to evaluate 
the density of states of each small string rigorously (with one end fixed and one end 
free) , 

Ui(£i) = e fe£ \ (2.20) 

In particular we find the same fin for all the strings. The total density of states is 
then given by 

j J] fadei) 5(e - J» = e toe , (2.21) 



we 



as required. 

It is now clear how to find the region occupied by the string. We measure this by 
determining the gyration, defined as the average size of the fluctuations of the free 
end from the fixed end, measured along null geodesies passing through the latter. 
Thermodynamics in the local inertial frames indicates that each small string has a 
spherically symmetric gyration with a radius ^Jel. To find the gyration of the large 
string we must (since they are average fluctuations) add those of the subsystems in 



11 



quadrature. Hence, regardless of which geodesic we choose, the combined system has 
a gyration VEi £i = yfe- 

Proceeding now to the volume dependence of u(e), we first make the usual quasi- 
equilibrium approximation that equilibrium is established much more quickly than 
any change in the metric so that the metric may be taken to be approximately 
constant in time when evaluating properties such as density. In order to simplify 
matters, we also make the (by now) familiar assumption that the metric is factorizable 
into parallel dimensions x and transverse ones y 

ds 2 = -n 2 dt 2 + g\\ ij dx i dx j + g± nm dy n dy m , (2.22) 

with the brane lying at y n = 0. We define parallel and total volumes as 

V\\{y') = j dxdy5{y-y')^/g[\ 
W\\{y') = J dxdy5{y-y')^g[\r](y,x) 
V = J dxdy^fg[\^fgZ 

W = J dxdy^gj ly /glr](y,x) (2.23) 

where 77 is a function that is one in the region of the random walk and zero elsewhere, 
and where g\\ and g± are the determinants of the metric in the brane dimensions and 
transverse dimensions respectively. Note that we are using x and y as shorthand 
for all the parallel and transverse coordinates. We also define an averaging over the 
extra dimension with an overbar, 



O = j^v^gfo) . (2.24) 
/ dy'Jgl 



For notational convenience we treat \Jg±(y) as a 'y-dependent transverse volume', 

V ± (y) = V&y) / dy', (2.25) 
so that the actual transverse volume is written 

VL = J dyJgAy). (2.26) 

Note that V = V\\V^_. From the discussion above, we have W\\ = e d wl 2 reflecting the 
fact that the string has typical size yfe. For example if we take a slice at y — y' this 
should be 

W\\ = e d ^ 2 = J dxdy5(y - y')^g[\r, = J dx^y'My') (2.27) 
so that i] must compensate for ^fg~\\ in order to make this volume independent of y'. 
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In order to examine the thermodynamics we now need to decide when we can use 
the microcanonical results. We first divide the parallel dimensions into small locally 
flat patches with ^/g± approximately independent of x. In each patch it is consistent 
to use the microcanonical results, provided there are no long range fluctuations. Thus 
we are generally prohibited from examining the NL systems since adjacent patches 
will not be in equilibrium. 6 

The random walk argument proceeds as before with one change. The factor that 
corrects for the translation of the walk in the excluded volume is not Vfi/VKy but 
rather V/W for the L[-l] system. In a toroidal compactification these are of course 
the same for the L[-l] case. Here however, we must take account of the fact that the 
walk is squeezed by the metric so that it may have more 'room' at one side of the 
compactified dimension than the other. This factor can also be written as 

V V\\ 

w = 4- (2 - 28) 

and we now find, for example, that the density of states in the L[-l] system is 

w(e) ~ (2.29) 
v± 

i.e. the same expression as in the flat space case but with all volumes averaged over 
transverse dimensions as in eqs.( |2.24| ) and ( |2.25| ). 



It is now possible to find log Z(j3, Vjj, V±) from uj(e) by integrating over e with a 
Boltzmann weighting, which we do for the various sytems in the following section. 
Note that the j3 appearing in the resulting partition function is the random walk 
definition of 'inverse temperature' whose precise physical interpretation is given by 

eq. (gjj . 



3. Stress-energy tensor in a bulk Hagedorn phase 

We now use these thermodynamic results to find the bulk energy momentum tensor 
during the Hagedorn regime. We may find the energy momentum tensor from 

*9 P " 



6 Also note that by examining the thermodynamics in a small region we do not artificially quench 
Kaluza-Klein modes as long as the local patches are much larger than the string scale. This is not 
the case in the perpendicular directions however because of the presence of macroscopic winding 
modes. Hence in a local patch with average volume Vji we may define an energy density by for 
example p = E/Vii where E is the total energy in the patch. However, for the thermodynamics, 
we must always include the whole transverse volume V± to avoid artifically unquenching winding 
modes. 



13 



where 0, V h V± are given by eqs. (tTTTD , (ggg) , (ggg) , (ggg) . We will treat the func- 



tional derivative with respect to g^ u in the following way. We assume that small 
changes in the metric correspond to making small changes in the volumes in logZ, 
e.g. for a single extra dimension 

¥-=l*>J z , f^- ^ 

ig 55 J SV L (x') Sg y , ( ' 

Then, in the case of only one extra dimension, we can write 

dyVgTs = -j^ = ^J d 5 x^r 5 , (3.3) 

and 

^ = ^— . (3.4) 

^55 2 v /^w/3 

Then from eqs. (|3.2| - |3.4| ) we can determine the functional derivative in eq. ( |3.1|) . Our 
ansatz automatically means that T 05 = and hence G05 = 0; in other words we are 
not considering energy exchange between the brane and the bulk. (In general there 
might be energy flux between the two.) 

We begin by evaluating the energy momentum tensor for our 'standard' case 
of L[— 1] (windings in all transverse dimensions) appropriate to high energies and 
densities; the resulting T^ v is presented in eq. flOll) . We then proceed to Lpy ^ -1] 
cases (in which windings are quenched in some of the transverse dimensions) and 
present the resulting T^ v in eq. (|3.27|) . 

3.1 Tp V at high energies and densities; L[-l] 



We now apply eq.(3.1) to the very high energy regime in which there are modes 
stretching across the whole space in all transverse dimensions. In the nomenclature 
of ref. these are the limiting L[-l] systems since d Q = 0. 

We first need the principal cosmological ingredient, the partition function. The 
expression for the density of states in the L[-l] system is 

u(e) ~ iefe", (3.5) 
v± 

which, when integrated with a Boltzmann weighting, gives 

log^, V„ VJ = + 0cM| _ V^-^M, (3.6) 

where a c and p c are a critical pressure and energy density (defined with reference to 
the brane dimensions, i.e. with dimensions E c /V\\) which are of order unity in string 
units 0. Here, subscript-c refers to critical quantities to remain in the Hagedorn 
phase. These are the successive terms in a saddle point approximation. 



14 



After doing the functional derivative we can replace (3 by using the saddle point 
result for the L[-l] system, 



^Ph(p-Pc), (3.7) 



V ± (f3 2 -f3 2 H ) ~\V X 

where as before, p = E/Vh. For the saddle point approximation to be valid we 
require || 

(p-p c )>i?i. (3.8) 

This is also the condition for the leading (first) term in log Z to dominate over the 
a c term, and for the a c term to dominate over the p c term. We neglect higher order 
terms. 

In order to write the resulting energy-momentum tensor, we define the parame- 
ters 

Vm V± 

a\\ = — ; ® ± = — (3.9) 

and a pressure 

P-i = ^L- (3.10) 

vr 

In the following section we shall argue that an « 1 wherever the Hagedorn phase is 
the dominant phase. If the off-diagonal components of the metric are small, then the 
variation with respect to g^ v gives the bulk (T^) terms; 

_ _~ P<x\\<*± 

T™ = -p_ x (2a|| - 1) « (3.11) 

where % labels the parallel x directions and m labels the transverse y directions. We 
have dropped the () notation for the expectation value of the energy momentum 
tensor, but shall continue to use hats to signify bulk properties such as the bulk 
density above, p. 

It is instructive to consider the conservation equations derived from the Bianchi 
identity Tf^ = 0. For example, restricting ourselves to the often discussed scenario 
with one extra dimension, y, and metric of the form 

ds 2 = -n 2 dt 2 + a 2 dx 2 + b 2 dy 2 , (3.12) 

the E-M conservation equations may then be written (we are using the same notation 
as in ref.0) 

^ + ( f Q -^ + (f °-f 5 5 )^0 

^- + it (- + -)+ -T?-3^7? = 0. (3.13) 
dy \n a I n a 
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We shall see shortly that a\\,a± ~ 1 so the first equation is a conservation law for 
entropy; 

& (l + y/Kh) - 2^W4 t (V,JjVZ) (3.14) 

where 

a = p + 2^/VZ. (3.15) 
Recall from ref.H that the saddle point dominance condition is 



^x/p«-^«l. (3.16) 

Hence to order 1/R"± we simply recover the expected entropy conservation law, 
- = 0. Note that the entropy is evenly distributed because of the dominant 
first term; indeed to first order Tq = a(y) where a is the local entropy density in the 
bulk. 

3.2 for the L[j ^ -1] systems. 

As the energy density of an L[-l] system falls below certain energy thresholds, the 
modes that are sensitive to the size of the transverse dimensions (e.g. winding modes 
in toroidal compactifictions) become quenched. Once this happens, the compactness 
of the quenched directions is of no further consequence for the thermodynamics, 
which aquires a different critical exponent, 7, and can even temporarily become NL 
(although the name 'non-limiting' is probably misleading for the reasons discussed in 
the Introduction and in ref.0). If there are many transverse dimensions, and they 
are anisotropic, there may be a few energy thresholds and 7 may assume several 
(increasing) values before the density drops below the Hagedorn density p c and the 
system finally leaves the Hagedorn phase and enters the Yang-Mills phase. 

In order to discuss the thermodynamics of these intermediate systems, we first 
introduce the function T]'(y) which gives support in the regions occupied by the 
random walking strings in the various transverse dimensions. (Note that rj(x,y) 
described the shape of a single string whereas rj'(y) is for all strings, so r)'(y) = 
indicates that strings do not extend this far into the bulk.) 

We define an averaging over the region given support by 7/; 



„ / dyrj y/gIO(y) 
I dy'rj'Jgl 



The single string density of states density of states now involves the ratio of the 
volumes 



V v , = J dxdyr]' ^/gZ^fg[\ 

W = J dxdyrfr}y/gly/g\\ (3.18) 
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(where rj is defined below eq. Q2.23Q ) so that 



Vrf _ V \W 

W W\\ ' 

The expression for the single string density of states becomes 

LO\£) ~ = L 

1 ; v ±v , ^ 



(3.19) 



(3.20) 



where in accord with our previous definition of r], we have assumed that the strings 
occupy a proper volume e d °/ 2 in the do dimensions where the modes do not fill the 
whole of space. As before we have j — d /2 — 1. 
Defining 

_ Vm 



'/' 



V ±rf 



(3.21) 



we may now use the expressions for the singular part of log Z found in ref . [ 

r(- 7 )7(/?-&) 7 , 7 ^z+u{o} 



log 3 



sing 



(-1) 



±gf f{P - ffl log(/3 - (3 C ) 7 € Z+ U {0}. 



(3.22) 



In the micro canonical ensemble, these expressions gives us the same relations as in 
eq. (|2.6| ) with p replaced by 

E 

7J==^, (3.23) 

II 7?' 

whence, substituting the appropriate expression for p from eq. fl2.6|) for each 7, we 
find the following for the L[7 7^ —1] systems; 



T ° = -V'P 



rpi. 



where 



p 7 



1h (2«n - 1 



=U log p 



' —Vjf> 1 windings 

no windings, 

if 7 = 4.5 



if 7 = 



(3.24) 



1, e -p 



Because, as we see later, ai 



if 7 = 1. 

1, we shall henceforth drop the overbar notation. 



17 



3.3 NL[7], WL and ML systems. 

For the NL systems, 7 > 1 and the microcanonical calculation tells us that 

(1 + 7) 

^-L-|. (3.25) 

In these cases we have already argued that they cannot be in equilibrium with the 
remaining systems since the temperature that we formally derive from the micro- 
canonical ensemble has T > T#. Indeed eq.( [2.11| ) also indicates the average total 
length of strings in the ensemble cannot be positive. The fact that there is no equi- 
librium solution to the Boltzmann equations is merely a different way of seeing that 
these systems are transient. Little more information can be gained from the ran- 
dom walk picture in this case. However it seems likely that, as in the flat space 
microcanonical calculation, these systems tend to lose energy to closed strings in a 
cosmological setting. 

The WL systems correspond to open strings that do not satisfy the saddle point 
condition. This is a cross over region where winding modes are just beginning to be 
excited, and the corresponding specific heat is therefore small (hence 'weak limiting'). 
Unfortunately this case is also difficult to analyse because the microcanonical ensem- 
ble does not agree with the canonical result (indeed there is no canonical equivalent 
of these systems). Consequently it is difficult to find anything that we may interpret 
as (IT) for these cases. The same is true for the ML systems which correspond to 
closed strings in the Hagedorn regime. Fortunately all of these systems have lower 
entropy than the limiting systems and can be neglected in the cosmology. 



3.4 Summary and discussion of T, 



We now summarize the results for the energy-momentum tensor. We first drop the 
overline notation of the previous section and simply redefine V± and V\\ to be the 
transverse and parallel volumes covariantly averaged over the region of the transverse 
dimensions covered by the strings. (See eqs. (|2.23| ),( |2.25| ),( |3.17| ).) We define an energy 
density p of strings 

p = V" (3 ' 26) 

For the limiting systems with critical exponent 7, L[7], we find that the 'bulk' com- 
ponents of the energy momentum tensor are given by 

1 o — P 

r p m ^ I — py transverse with windings (3 27) 

" 1 transverse without windings, 
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where p = p/V± and 



i 



if ■y = —I —I I 
11 f ±, 2 , 2 



V", 



3/2 



logp 



if 7 = 



(3.28) 




-p 



if 7 = 1, 



± 



wherever there are strings present, and zero otherwise. The approximation in the 
transverse components is valid when the parallel volume changes only by a small 
fraction over the transverse directions. 

For the remaining NL, WL, ML systems we do not know how to consider the 
cosmology or indeed whether they have any meaning at all in a cosmological setting. 
In the NL case large scale fTuctations make it impossible to have a stable macroscopic 
system. For the WL and ML systems, one cannot properly take into account the 
effect of winding modes. The fact that the canonical and microcanonical ensembles 
disagree in a toroidally compactified space tells us that they play a dominant role in 
the thermodynamics, but it is not possible to formulate the required microcanonical 
treatment in cosmology since well defined global parameters such as total energy are 
generally lacking. 

As expected T ° resembles the local energy density of strings. The T\ represents a 
relatively small pressure coming from Kaluza-Klein modes in the Neumann directions 
and T™ is a negative pressure coming from winding modes in the Dirichlet directions. 
If we T-dualize the Dirichlet directions these 'winding modes' also become Kaluza- 
Klein modes in Neumann-directions and T™ becomes positive. Thus negative 
reflects the fact that we have T-dualized a dimension much smaller than the string 
scale thereby reversing the pressure. For this reason negative Tjf is expected to be a 
general feature of space-filling excitations in transverse dimensions. 

Those who are familiar with T-duality may suspect an apparent conflict with 
this reversal of pressure when we T-dualize the extra dimensions. However we stress 
that T-duality is maintained. Pressure is merely defined as the change in free energy 
with increase in volume. Thus since free energy is invariant under T-duality the 
pressure must reverse sign. On the other hand the cosmological consequences in the 
T-dual system must of course be the same. In particular gravitons do not propagate 
in the extra dimensions of the T-dualized system (i.e. their Kaluza-Klein modes are 
heavy) since the extra dimensions are much smaller than the string scale. However 
gravitational degrees of freedom come from the closed string sector in the bulk which 
has both Kaluza-Klein modes and winding modes in all directions. Hence in the 
T-dual system the light gravitational degrees of freedom in the extra dimensions are 
winding modes. The original transverse components of Einstein's equations (which 
are valid only in the approximation that the Dirichlet directions are much larger 
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than the string scale) must in the T-dual system be replaced by the equations of 
motion derived from the action now expressed as a sum over gravitational winding 
sectors. The upshot is that T-duality dictates that the resulting equations have the 
same consequences for the brane cosmology as the original system. 

4. Cosmological Equations in D = 5 

In the previous sections we obtained the bulk stress-energy tensor at temperatures 
close to the Hagedorn temperature for the most general limiting cases in arbitrary 
numbers of dimensions. We now proceed to a discussion of the cosmology for which 
we will for the most part consider a 'toy-model' system with the following restrictions. 

The first restriction is that we implicitely be considering a D-brane configuration 
that has 3 large parallel dimensions (i.e. the 'observable universe') and only one 
transverse dimension y that supports winding modes. Such a configuration could 
possibly arise in a 5 dimensional intersection of D-branes. A more realistic case would 
probably include all D = 10 dimensions but the Einstein equations are significantly 
more complicated in this case and we leave their full examination to a later paper. 
However in section 5.6 we will discuss why we expect to find the same behaviour in 
higher dimensions. 

The second restriction we make is simply for definiteness; we assume adiabaticity 
when solving Einstein equations. We shall find in this section and in section 5 
both power law and exponential inflation. In principal, as we will see in section 
6, it is therefore possible to solve the horizon problem even with adiabaticity. In 
practice however it is not possible to introduce a priori the entire entropy of the 
observable universe (i.e. S ~ 10 88 ) onto a primordial brane of order the string scale 
without either destroying the brane or having an unreasonably small string coupling 
(i.e. g s ~ 10 -88 !). In section 6 therefore we shall discuss possible deviations from 
adiabaticity. 

We use the metric of eq.( |3.12|) , 

ds 2 = -n 2 dt 2 + a 2 dx 2 + b 2 dy 2 , (4.1) 

which foliates the space into flat, homogeneous, and isotropic spatial 3-planes. Here 
x = xi,X2,x% are the coordinates on the spatial 3-planes while y is the coordinate 
of the extra dimension. For simplicity, we make a further restriction by imposing Z 2 
symmetry under y — > —y. Without any loss of generality we choose the 3-brane of the 
'observable universe' to be fixed at y — 0. (Even if the brane is moving we can always 
reparameterize the theory so that the 3-brane is fixed in coordinate space, although 
in doing so we lose further freedom of gauge choice []I2"| .) With our assumption of 
homegeneity we can associate a scale factor with the parallel dimensions a(t, y) and 
one for the 'extra' transverse dimensions b(t,y). We define ao(t(r)) = a(t,0) as the 
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scale factor describing the expansion of the 3-brane where t(r) = / drn(T,y = 0) is 
the proper time of a comoving observer. 

Several authors |T(J - jOJ have presented the bulk Einstein equations but for 
completeness we briefly restate the results; 
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where we use the notation of ref. 10 in which 



k 2 = 87rG = 8tt/M 5 3 , 



(4.3) 
(4.4) 
(4.5) 

(4.6) 



where M 5 is the five-dimensional Planck mass and the dots and primes denote dif- 
ferentiation with respect to t and y, respectively. Note that, as stated earlier, our 
ansatz implies that T05 = in the bulk. 

The bulk equations only apply in an open region that does not include the 



boundary. It is the Israel conditions which connect the boundary and the bulk ]10| 

t 



J2 (K^-Kh^) 



"IAV 1 



(4.7) 



ifaces 
curvat 

boundary surface and we have defined 



where is the extrinsic curvature, the sum over faces is for each side of the 



t 



2 6 S 'boundary 



(4- 



as the energy momentum tensor on the boundary. We will assume that this energy 
momentum tensor on the boundary can be written in a perfect fluid form. For 
example, 



^0 — Pbr 



and 



-Pbr 



(4.9) 
(4.10) 
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where /v and are the energy density and pressure, respectively, measured by a 
comoving observer. 



We can now find the Israel conditions specific to the metric in eq. ( 3.12 ) . Because 
we are considering a brane at a Z 2 symmetry fixed plane, from the sum over faces we 
get two identical terms, merely resulting in a factor of two. In a more general space- 
time, the Israel conditions would still hold, but one would have to explicitly add the 
contributions from the two sides of the boundary surfaces since these contributions 
would no longer be equal. With the Z 2 symmetry, the Israel condition becomes, 

3[a'/a] = -k 2 b p br 

3[n'/n] = k%(2 Pbr + 3 Pbr ) (4.11) 

The same result can be found by adding 5(y)t^ to Tjf and equating the coefficients 
of 5(y) in Einstein's equations. 

It is interesting at this point to consider the relative contributions of the stringy 
excitations and of the D-brane tension to the cosmology. In particular we, somewhat 
counterintuitively, find that the diffuse stringy component can have a dominant effect 
on the cosmology even in the weak coupling limit. To see this, first note that the 
effective gravitational coupling k is given in terms of the string coupling g s by 



1 m?-*Vol D _ 5 3 

— M 5 . (4.12) 

K 9 S 



where VoIb-^ is the volume of the compactification from D = 10 to 5 dimensions, 
m s is the string scale and M 5 is the 5-dimensional Planck mass. The intrinsic tension 
of a D-brane is given by 

T Dp = — \— ~ - * (4.13) 

{2ir)Pg s k^VoId-s 

For convenience we will henceforth assume that VoId-5 ~ m^ 5 so that g s and k are 
of the same order of magnitude in string units (where m s ~ 1). (This need not be 
the case if some of the other space dimensions are compactified with a radius much 
larger or smaller than the string scale.) Thus the intrinsic tension of the D-brane 
satisfies p br ~ and in the small coupling limit one might expect it to dominate 
the cosmology. However, this is not the case. If we substitute the Israel matching 
conditions into the cosmological equations ( |4.2| - |4.5| ), we see that the contributions 



of the D-brane tension and of the stringy components in the cosmological equations 
are of order 

-4 2 i -2 P -2 - 

k p br and k — — , k p 7 
V± 

respectively. When the density of the string gas is close to the critical Hagedorn den- 
sity, p c ~ 1 (the effective lower bound) , all the terms are of order k 2 and are therefore 
comparable. Since the density of the string gas can be significantly larger than the 
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lower bound of p c ~ 1, clearly the bulk stringy contribution can be significantly 
larger than the brane contribution. 

We can further see the dominance of the bulk stringy component in the weak 
coupling limit. All our work has assumed that 

P^T Dp . (4.14) 

If this condition is violated, the thermal energy of the brane is larger than its rest 
mass and one would expect our perturbative treatment of the D-brane to break down 
(see ref.|| for possible outcomes.) (Note that this bound does not apply directly 
to orientifolds which are stuck at fixed points.) Since Td p oc -, we can see that 
the amount of thermal energy that can be loaded onto the brane increases as the 
coupling gets weaker. When this bound is saturated, the stringy contribution to the 
cosmological equations can be O(k) while that from the intrinsic tension is 0(k 2 ); 
thus weak coupling is advantageous for string dominance in the cosmology. This 
rather surprising fact will become clear in the examples of the next section. 

5. Results: Behaviour of scale factors a(t) and b(t) 

In this section we present results for the behaviour of the scale factors a and b. As 
discussed above, we will for simplicity and definiteness mainly consider 5 dimensional, 
adiabatic systems. In sections 5.6 and 6 we will consider how it may be possible to 
relax these restrictions. 

Our results are obtained by solving Einstein's equations together with the con- 
straints provided by the Israel conditions. We use the energy momentum tensor 
derived above appropriate to a primordial Hagedorn epoch. First we will consider 
the regime Lpy] under the assumption that the brane is energy /pressureless and that 
5 = 0. We also discuss some constraints imposed by the stability of the extra dimen- 
sion. We then catalogue more general types of behaviour that can occur for b ^ 
and discuss a semi-realistic example in which b is constrained. 

5.1 Hagedorn inflation; A^ r = ^buik = b = 

Let us first impose 6 = simply as an external condition; i.e. the extra dimension is 
fixed in time. We also set the cosmological constants both on the brane and in the 
bulk to zero, Af, r = A^^ = 0. All other contributions to the energy momentum local- 
ized on the brane (for example massless Yang-Mills degrees of freedom or possibly an 
NL subsystem) have an entropy that is subdominant to the limiting bulk degrees of 
freedom as long as p > p c (where p c is a critical density of order 1). For the moment 
we will therefore neglect them and set a'(t, 0) = n'(t,0) = 0. In this discussion the 
brane at y — is playing no role in determining the evolution of the cosmology; the 
scale factor a changes purely as a result of the bulk equations. (We shall see shortly 
that consistency of the full solutions requires a second brane at y = I). 
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A. Solutions to the T 55 equation 



It is simple to solve the 55 equation for the scale factor a(t, y) at the origin, 

a(t,0) = a (t), (5.1) 

for all the L[y] systems. 

Al. L[-l]: In the L[-l] systems (the high energy case with windings in all transverse 
dimensions) we see that T| ~ y/p ~ a -3 / 2 (where, again, p is the local energy density 
measured with respect to the volume Vm) and hence, 



a 



(t) ~ t 4 / 3 , (5.2) 



i.e. power law inflation. 

A2. L[0]: We find that 7 = gives us a period of exponential inflation. To find this 
solution we write a 2 , = exp F and the T| equation becomes 

F + F 2 = ^(logp(0)-^F 
which may easily be solved to give 



k? 



j-log(p(0)e3/*) (5.3) 



where p(0) is the initial density at t — 0. Initially the second term in the exponent 
dominates and there is exponential inflation. This solution has an automatic end to 
inflation when dao(t)/dt ~ 0, i.e. when the two terms in the exponent are roughly 



of the same magnitude. Hagedorn inflation ends at time t ~ k 1 vV_l logp(0) which 



corresponds to log ^| ~ logp(0). During a period of adiabaticity, the entropy 
and energy density dilute as 1/a 3 , because the Hagedorn phase is always to a first 
approximation like pressureless matter (p_ x <^ p). Hence we find that the above 
condition for inflation to end happens when p(t) ~ 1 (in string units), just as the 
system is dropping out of the Hagedorn phase. 

A3. Other systems: We summarize the energy momentum tensors and cosmo logical 
behavior of ao(t) for the remaining systems in table 1, where p = 3 in our 'toy- model'. 
In all these cases it should be noted that during a period of adiabaticity, the entropy 
density dilutes as 1/a 3 ,. The reason for the various different types of behaviour is of 
course due to T| which may drop off more slowly than 1/a 3 , and in the L[0] case drops 
only logarithmically with the expansion. Note inflationary (superluminal) expansion 
for L[ 7 < 0]. 

In adiabatic systems the Hagedorn regime and hence the inflationary behaviour 
we have found eventually come to an end. For a system to be in the Hagedorn 
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regime requires an entropy density higher than the critical Hagedorn density (of 
order 1 in string units). Below this density the energy momentum tensor and hence 
the cosmology is governed by the massless relativistic Yang-Mills gas (the gas is 
present on the brane even in the Hagedorn regime but is sub dominant). Thus there 
is no problem exiting from the inflationary behaviour. In fact the main issue is how 
long inflation can last. We return to this question later when we discuss how inflation 
can be sustained. 
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Table 1: Cosmological regimes for open strings in the Hagedorn phase with p large parallel 
dimensions. The constants C and D are given for a p = 3 (i.e. the 'observable universe' as 
a 3-brane) in eq.([5.3|). Note inflationary expansion for L[7 < 0]. 



B. Solutions to the Too an d Ta equations 

We can get an approximate solution for the remaining equations which is valid 
in all the regions of interest as follows. For small y we write 

y 2 

a(t,y) = ao(t) + —a 2 {t) + ... 
V 2 

n(t,y) = l + |-n 2 (t) + ... (5.4) 

where henceforth we normalize n(t, 0) = 1. First, given our constraints above, the 05 
equation is trivially satisfied. As before we solve the 55 equation with p 7 = V^ 3 ^ 2 p 2 ^ 3q 
and find 

a (t) w At q , (5.5) 

where 

h y A 2/* = * a ( 5. 6 ) 

° 3g(2g-l) 1 J 

Solving the 00 and ii equations with Tq = —k 2 p/V± = —k 2 p we find 

a 2 (t) = bla 



, 2 (<£_ tr£ 

t 2 3 , 



Mt) ^ b i(^ + ^i±m. (5.7) 



t 2 
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Phase diagram constraints on the transverse dimension and an ~ 1 



One important aspect of the approximate solution is that it is valid in all Hage- 
dorn dominated regions. Obviously the approximation requires that the y 2 terms are 
small. Looking at a(t, y) we see that the condition is a 2 y 2 < a or 

k 2 p < JL (5.8) 

Recall from ref. @] that the Hagedorn regime occupies a small-coupling region between 
the Yang Mills phase and the black brane phase. Thus if the energy density is too 
high, thermodynamics is dominated by black branes and our entire perturbative 
derivation ceases to be valid. In addition there is a 'holographic' constraint which 
is saturated when the entropy is equal to that of a black Dp-brane filling the entire 
transverse volume (for a review see ref. [ 



k 2 pV ± « k 2 p < Vf~*-*, (5.9) 



where as before (see eq.( |4.13| )) we assume for convenience that the volume of the 



compactification to D dimensions is ~ 1 so that k ~ g s . In the case that the 9 — p 
transverse dimensions are isotropic, this is the same as eq.( |5.8|) . The present case is 
(for convenience) taken to be one in which there is only one transverse dimension of 
size R± which is much larger than the string scale. We then set D = 5 and p = 3 
and again find eq.( |5.8|) . Thus our approximate solution for the particular case above 
breaks down precisely where the holographic bound is saturated - when the event 
horizon of a black-brane fills the transverse dimensions. 

Note that, depending on the dimensionalities, the entire calculation (i.e. based 
on Hagedorn regime thermodynamics) can already break down at much lower den- 
sities corresponding merely to black brane dominance. Considering for defmiteness 
black-hole type manifolds with planar asympotics, Si x T£, it was seen in ref.0 that 
when 

/t 2 p>l (5.10) 

the dominant component in the thermodynamics (i.e. that phase with the largest 
entropy) is the horizon of a Schwarzschild-like black brane. Since R± is larger than 
the string scale, eq.( [5.10D can be satisfied at a density much lower than the density 



of the holographic bound above. 

Adding in the constraint that the Yang-Mills entropy be less than the Hagedorn, 
we find that the phase diagram restrictions can be summarized as 



mm 



in{^- 3 - p ,l} > £ 2 p> k 2 N 2 (5.11) 



where N is the number of branes situated at y — 0. (The additional requirement of 
P ^ Pc ~ 1 explains why we need to be at a small string coupling, k 2 1.) 
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One consequence of the phase diagram constraints is that in the full solution a 
is not allowed to vary significantly in the range y G [—£,?], and hence, as promised, 
Q!||(0) ~ 1 in all the regions of validity. Note that a" and n" must still compensate 
for Tq being greater than T| (the saddle-point dominance condition) and these two 
requirements are compatible only at weak coupling. 



Stabilization and the second brane 



A nice check of our approximate solution is that it should obey the stabilization 
condition discussed in |10|, [L4| for the case where T 05 = (as we are assuming 
throughout), 

/ dy^{T° + T\ - 2T 5 5 ) = 0. (5.12) 

To see that it indeed does without any further assumptions we need to consider 
what happens to the solution at the boundaries of the compact space. Assume that 
the large transverse dimension is compactified on a range [—1,1] and recall that for 
convenience we are assuming Z 2 symmetry under y — > —y. Thus at y = I the 
solutions have a discontinuity in the derivatives of a and n which is effectively a 
second brane. Using the discontinuity condition at y — I, 

k 2 

a'_i - a' +l = — —Pbra(l)b 
k 2 

- n' +l = +y(3p 6r + 2p br )n(l)b , (5.13) 

we find that this corresponds to additional terms in the energy-momentum tensor 
given support at y = I by a delta function; 

AT ° = 5(y - l) Pbr = S(y - l)R ± (p - 3- T ^ 



k 2 t 2 J 

at* = <% - i)p» = Kv - QR± (-P7 - ' (5 - 14) 

where R± = b l and where we have neglected terms of order k A p 2 R\ in accord 
with our previous discussion. Note that the signs of the brane contributions are the 
opposite of those coming from the bulk, so that eq. (|5.12| ) now reduces to 

dyft^.+T*\=0. (5.15) 

This equation is automatically satisfied since it is the 55 component of Einstein's 
equations. 

In the stabilized system (i.e. with the 6 = ansatz) the brane we find at y = I 
inevitably has a peculiar equation of state that does not resemble a Yang-Mills or 
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cosmological constant on the brane. Nevertheless the stabilization condition should 
always be satisfied so that, in a full model which includes a stabilization mechanism 
(given our assumptions for a' = h' = n' = at y = and Z2 symmetry), the 
gravitational sector should behave as if there were a brane with these characteristics 
located at y — I. 

5.2 L[ 7 ] with %k\ T - A bulk = and b + 

We now study the case where the extra dimension is not stabilized but there is still 
no nett cosmological constant. By 'nett' we mean that the contribution from a 1 and 
n 1 in G| cancels the bulk contribution on the RHS of the Tj equation. The condition 
for this is 

k 2 

Kett = -JF, A lr + A bulk = 0. (5.16) 

In this case we could again, by studying the 00 and ii equations, try to find an 
approximate solution which is valid throughout the Hagedorn regime as we did in 
the previous subsection. However we would not learn anything further by doing this 
since we could always independently adjust a 2 (t) and n 2 {t) to satisfy the equations. 
(There is nothing other than the phase diagram constraints to predetermine these 
parameters.) In addition eq. ( |5.12| ) need not be satisfied (and in fact the integral is oc 
b), but one should keep in mind the fact that the full solution again has discontinuities 
in a' and n' at y = I so that in general there should be a brane (or something that 
can compensate in a similar way) situated at y = I even when b 7^ 0. The next 
subsection shows that by eliminating the second brane at y — I one gains an extra 
constraint on a" or n" and hence useful additional information from the ii and 00 
equations. 

Therefore we will first focus on the 55 equation. Since a ~ 1, for 7 = —1, — |, ^, 
the pressure is 

-37 ^3 

Pj oc aT- 1 b 2 . (5-17) 

We will now catalogue some possible types of behaviour for a with different ansatze 
for b. 

Taking A nett = 0, we find a family of power law solutions for 7 = — 1, — |, | of 
the form 



a (t) w At q 

6 (t) w Bf (5.18) 



where subscript-0 indicates values at y = 0, 



7-14 

37 3 k 2 , 
A~Bi = — -, 5.19 
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and where r is arbitrary. Note that if we choose r = 4/3 then b ~ t 4 / 3 and a is 
stationary whatever the value of 7 (the complementary situation to b stationary and 
a undergoing power law inflation in eq. (|5.2|) ). 

We also find a family of hyperbolic solutions of the form 



a (t) /sinh2C(t + t 



ao(0) V sinh2C£i 
b {t) 



/sinh2C(£+£i) \ 1-7 :r _ _i _1 1 

V sinh2Cti J I ' 2> 2 



6o(0) 



log f^2C( Wl )\ if q 

iU & V sinh2Cti ^ 



,-p(S|fi^) if 7=1- 
where C and ti are defined by 

C 2 = ^> 7 (0) 
a (0) 



6 

Ccoth2Cii, (5.20) 



a (0) 

p 7 is given by eq. (|3.28| ), and where e.g. p 7 (0), p(0) = p 7 (t = 0), pit = 0). Note that 



the solutions with ti — > ±00 describe the scale factor exponentially increasing or 



decreasing; e.g. for 7 = —1, — |, ^ in this limit we have 



qp(t) 
a (0) 
60 (i) 



exp (±Ct) , ti — > ±00 

,„„„, - exp (±^Ct) , ti -> ±00. 



(5.21) 



In the hyperbolic solutions T 5 5 is constant in time (see eq. (|5.17| ). Both expansion 
and collapse of a are possible, since b compensates appropriately (in such a way as to 
keep T| constant). It seems rather paradoxical that the compensation can go either 
way depending on the sign of 7. For example, in the 7 = 1/2 solution of eq.( |5.21|) , if 
a(t) is expanding then, in order to keep T| constant, b has to expand exponentially 
as well. The solution in the previous subsection for the b = 0, 7 = system (section 
5.1.A2, eq.(p3|)) is similar, although there is a slight difference because there b is 
constant, resulting in the e~* 2 'exit' term. We should also remember (bearing in 
mind our previous discussion) that if adiabaticity is assumed then the density T° is 
always diluted exponentially fast as T° oc exp (— Ct~^j . Finally we should add that 
which, if any, of these solutions is appropriate depends on the initial value we choose 
for b . 

5.3 L[ 7 ] with ^A 2 br - A bulk ± and b ^ 

Next we consider the case of additional cosmological constants in the brane and bulk. 
We transfer the a' and n' terms to the RHS of the 55 equation and write 

T* eff = -(p, + A nett ), (5.22) 
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where 



An 



ett 



A. 



bulk 



12 



-A 



br • 



(5.23) 



First we can see (rather trivially) that the hyperbolic solutions of the previous 
subsection in eq.( |5.20|) still exist if we modify C to C", where 



C 



12 



6 



(A nett + p 7 (0)) 



(5.24) 



and provided C' 2 > 0. Not surprisingly we recover the usual cosmo logical constant 
driven inflation when we set p = 0. When C' 2 < (which is the case for A nett < 
and |A nett | > p 7 (0)), we find a singular collapsing solution for a; 



ao(t) _ f sm2\C'\{t + t x ] 
a (0) " 

bo(t) 



6o(0) 



where t\ is defined by 



wx2\C% 

n2\C'\(t+t-_ 
sin2|C|*i 

, / sin2|C'|(Wi) \ 
10 & V sin2|C'|ti J 



^ sin2K "!■/-'-■' i i \ 



\ sin2|C'|ti 



if 7 = — 
if 7 = 
if 7 = 1. 



1 -± i 

' 2' 2 



flo(O) 
ao(0) 



IC'I cot2|C"|ti. 



(5.25) 



For the power law solutions eq. (|5.18| ), we note that the Hagedorn contribution 
to T| varies as 1/t 2 and therefore decreases in time compared to the constant A nett 
contributions. Thus we can consider the 55 equation in the two limits that T| e/ -, in 



eq.( 5.22j) is dominated either by the cosmological constant term or by the Hagedorn 



term. For example, if 7 



_i 1 

'' 2' 2' 



7:, we have 



Oo(t) « A(t-to) 9 
6 (t) &B(t-t ) r 

a (t) _ / sinh2C"(t+ti) \ 2 

a o(0) ~~ V sinh2C"ti / 

bp(t) _ / sinh2C'(t+ti) \ 7^7 

^0 (0) ~~ V sinh2C"ti ) 

a (t) _ / sin2|C"|(t+ti) \ 2 

Oo(0) ~~ V sin2|C'|*i j 
bp(t) / sin2|C , |(t+ti) ^T^7 
feo(0) V sin2|C'|ti 



2 5eff 



I I 
\ 1 5eff\ 



I T - ^ I 
\- L 5eff\ 



-Pi 



\Kett\ >\Pj\, C' 2 > 



|A neit | » lp 7 |, C' 2 < 



(5.26) 
(5.27) 

(5.28) 



where to is the time at which a is released from its small initial value. 
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5.4 Bouncing universe with softened singularity 



If A neU < the above solutions admit an oscillating universe. For example if r < 4/3, 
we find a universe in which a negative cosmological constant causes a to follow the 
sin curve of eq. (|5.28| ) which reaches a maximum at \C'\(t-\-t\) = tt/4 before heading 
towards zero. However, unlike the case of pure negative cosmological constant in 
ordinary FRW, we do not hit a singularity because, when the scale factor a becomes 
small, A nett + p 7 in eq. ( |5.22| ) becomes positive and T| is dominated by p 7 . Instead 
of hitting a singularity, a reaches some minimum value and rebounds upwards on 
the power law solution of eq.( |5.26| ) (with to now marking the time of the rebound, 
d(to,0) = 0). The scale factor a follows the power law curve until the Hagedorn 
contribution to T| (which varies as 1/ (t — to) 2 ) drops below the cosmological constant 
piece, at which point we pick up the sin curve of eq. ( |5.28|) and have completed a single 
oscillation. The oscillation is around the completely static solution (with a = a = 
b = 0) where the pressure p 7 is exactly compensated by a negative cosmological 
constant, 

p 7 = -A nett . (5.29) 

The behaviour of b depends on the value of 7. For 7 = | we find that b contracts 
and expands with a. In this case, at small a we must have a power law solution 
with r > 4/3 so that the contracting/expanding power law solution matches onto a 
contracting/expanding sin solution. On the other hand for 7 = —1 or — ~, b contracts 
and expands out of phase with a, and we require < r < | at small a. Clearly if, at 
small a, r « then from eq.( |5.26|) b is almost static. 

We stress that this singularity smoothing behaviour may or may not occur de- 
pending on the initial conditions (e.g. b) and additional constraints on the system. 
For example, it does occur in numerical solutions with 6 = 0. 

It is encouraging to see that one contribution of stringy physics to the cosmology 
can be to soften the singularity that would otherwise appear at a = 0. This is a 
familiar aspect of string theory, but an appealing feature of the Hagedorn phase is 
that we can find it purely perturbatively. 



5.5 Full solution with an example of a 'physical' brane 

In the previous subsections we discussed the cosmological behaviour for various 
ansatze. However, apart from the 6 = constraint (which can be motivated by 
some unknown stabilization mechanism) the other ansatze are unmotivated, and we 
have still to show that the behaviour can arise in 'realistic' physical systems. The 
reason for this unwanted freedom is that the system of equations is underconstrained 
- there are four parameters (a , b Q , a 2 , n 2 ) but only three independent equation 
(00, ii, 55). Therefore, we now consider an example in which we impose an addi- 
tional constraint coming from a particular choice of equation of state on the brane. 
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This allows us to find a full solution, in which the behaviour observed in the previous 
subsection arises in certain limits. 

As we saw when we discussed stabilization, there is a discontinuity in a' and n' 
which corresponds to a brane located at y = I. Generally, the equation of state for 
this brane will look rather unrealistic (as indeed it does for the static b = solution) 
and we are forced to account for it by invoking unknown contributions from the 
gravitational sector. Our first assumption is therefore that there is no second brane 
(or rather, in view of the need to conserve Ramond-Ramond charge, that only the 
brane at y = has any cosmological effect). 

With this assumption, continuity at y — I requires a\ = n[ = and we must drop 
the approximation that the single brane at y = (our would-be universe) is 'empty' 
(in the sense that a f — n' — 0) and consider its energy and momentum as well. The 
Israel conditions tell us that this energy density is given by 

3[a'/a] = -k 2 b pbr 

3[n'/n]o = k 2 b (2p br + 3p br ). (5.30) 

We now make an additional assumption about the equation of state of the brane at 
y = 0; we assume it is also Hagedorn-like, i.e. as well as the bulk L[y] degrees of 
freedom there is an additional localized Hagedorn component with excitations near 
the brane. Even though we do not understand the thermodynamic behaviour of 
this subsystem, we can be fairly confident that its energy momentum tensor is like 
pressureless matter (as all the Hagedorn systems are to a first approximation - they 
all have \pbr\ ^ \Pbr\)- Note that, to simplify this example, we will neglect the brane 
tension contribution which acts like h br . 

For p br pa it is therefore sufficient to impose the relation 

n'/n = -2a' /a (5.31) 

for all y. Integrated, this gives 

< 5 ' 32 > 

Together with a[ = n[ = this means that the approximate solutions for a(t, y) and 
n(t, y) are going to be of the form 

a(t, y) = a (t) + <n(t) (^\y\ - 

•*•»> =i - 2 S( w -9- (5 - 33) 

The 05 equation may easily be solved and we find 

b\(y) = a'a 2 , (5.34) 
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where X(y) is a constant of integration. Note that this and eq. ( |5.30| ) imply that 

--TPbr = ~77\3 ( 5 - 35 ) 
3 a (ty 

so that the brane density is indeed diluted by the expansion of the scale factor a like 
pressureless matter. 

We now apply this to the 00, ii and 55 Einstein equations which now are more 
constrained by the continuity at y — I and by our choice of equation of state. Some- 
what fortuitously the condition in eq.( 5.31| ) gives a cancellation of the a" and n" 
terms and the equations at y = reduce to 

oo/oo + &o\Ao_ _^o_ _ * 

oo \a>o b J «o a o b o 3 

3 A| + ao /do + M + 2 ao + bo = _ k2 

% a \a b J a b 

— + — +— = +^-P 7 - ( 5 - 36 ) 
ag \aoJ a 3 

Dominance of the p pressure term in the 55 equation is only possible if 

\ 2 o/al < (a/ a) 2 ~ \k 2 p 7 \ < K 2 p (5.37) 
and so the 00 equation must be dominated by the A' term. Solving gives 

_ i « i.» (5.38) 
Consistency of our solutions therefore requires that p& r £ p and hence by eq.( |5.37|) 

£ 2 p « & « _L (5.39) 

P Kj. 

As for the previous solutions, this corresponds to the holographic constraint in the 
phase diagram, and implies small coupling. Thus we again observe that, because 
only £ 4 Pbr appears in the Einstein equations, the cosmology can be dominated by 
the 'bulk' energy-density even when p < pt, r , provided that the coupling is small. 

For the L[-l] systems the ii and 55 equations may be solved by substituting 
d = ao , ( nK and c = so that 

a (0)6 (0) ao(Or 

c 2k 2 _ d Ak 2 _ Ak 2 



-P-i 



c 3 d 3 3 

Defining 



P-i = — — p-i(0)d~ 3/2 . (5.40) 
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we find power law and hyperbolic behaviour emerging in different limits. When 



3 

!i)2 



we have power law behaviour (eqn. pMS ) with 



* - *0 » TU /n^ 2 ( 5 - 41 ) 



So for large r ps 4/3, 5/6. On the other hand when 



3 

* - *o < (5.43) 

K 4 p_l(0) 2 



we have \fd ~ — w / k A p_i(0) 2 ps yd(0) and hence we find the hyperbolic solution of 
eq.( |5.27| ) with expanding a and collapsing 6. 

5.6 More dimensions 

Finally we briefly remark on the extension to higher numbers of dimensions. In type 
I,IIA/B models the total number of dimensions is D = 10 and the energy momentum 
tensor is therefore of the form = diag(TQ,p 7 ,p 7 . . . , — p y , — p 7 . . . , . . . , 0). A Dp- 
brane has 9 — p large Dirichlet directions giving 

-1 < 7 < min(l,^-^l . (5.44) 



2 

As we discussed in the introduction, it is always true that 7 < 1 because for p < 5, 
the only possible limiting system has excited winding modes and is the 7 = — 1 
system. 

A diagonal metric with universal scale factors is consistent with the above form 
of T". We choose the metric to be a function of a single transverse dimension 
which we label y, and add D — 2 — p additional transverse dimensions labelled y m , 
(m = p + 2, ..£>); 

ds 2 = -n 2 dt 2 + a(t, yfdx 2 + b(t, yfdy 2 + d(t, y) 2 dy m dy m . (5.45) 

This metric gives 4 independent equations (00, ii, 55, mm) for the higher dimensional 
systems. Most importantly, a", n" and d" do not appear in the 55 components of 
[|n]]. Consequently the 55 equation may again be treated separately from the 00, 
ii, mm equations. The latter 3 equations may always be solved independently (as in 
the D = 5 case studied in detail in this paper) using the three unknowns a", n" and 
d" . Thus, although we will not discuss the higher dimensional case in detail here, we 
can immediately see on dimensional grounds that the solutions to the 55 equations 
with d — b — must be as shown in table 1 (with the appropriate value of p) . An 
amusing aspect of this is that power law inflation in L[-l] systems [i.e. the universal 
high energy system) requires p < 3 regardless of the total number of dimensions. 
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6. Sustaining inflation and solving the Horizon Problem 



One of the most striking features of the cosmological solutions we have found is 
that they automatically predict a superluminal growth of the scale factor which is 
usually an ingredient of standard inflation. The remaining point that needs to be 
addressed is that our analysis assumed adiabaticity, and in order to realistically solve 
any cosmological problems, this assumption (as we will see) almost certainly requires 
modification. In this section we discuss solving the horizon problem beginning with 
the causality condition and then addressing the issue of adiabaticity. We then discuss 
how nonabiabatic effects may contribute to a sustained period of inflation. 

6.1 Causality Condition: 

An estimate of the size of the observable universe today is given by the distance light 
could travel between photon decoupling and now, 



Note that d b s = 0(1) x (t — td ec ) for a oc t p and p = 0(1) between tdec and t Q . Here 
a is the scale factor. We can compare the comoving size of the observable universe 
to the comoving size of a causally connected region at some earlier time t p : 



The observable universe today fits inside a causally connected region at t p if 



Here, subscript-o refers to today and subscript-p indicates some primordial time 
before the inflationary period of interest. If condition eq. (|6.3| ) is met, then the horizon 
size at t p (before nucleosynthesis) is large enough to allow for a causal explanation of 
the smoothness of the universe today. Note that more creative explanations of large 
scale smoothness may not involve comparing these two patches. For example, in the 
context of the brane scenarios, one might imagine that two regions of our observable 
universe which seem to be causally disconnected might in fact have talked to each 
other because of a geodesic between them that went off our brane, into the bulk, 
and then back onto our brane at some distant point, as demonstrated by Chung and 
Freese [18]]. In the remainder of our discussion here we restrict ourselves to the case 
where eq. ( |6.3| ) is relevant; this is certainly the case for all the brane and boundary 
inflation models considered to date. 

For power law expansion of the scale factor both before t p and after £<&, c (which 
may or may not be the case), we can take t ~ H~ l during these periods. The 






(6.2) 




(6.3) 
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causality condition eq.( |6.3| ) then becomes 

1 



r h ~ TnT' (6 - 4) 

or, equivalently, 

Q>p ^today 

We take the Hubble constant today to be 

H = al /2 T 2 /m pl (t ), (6.6) 

where a Q = (8ti / 3) (it 2 /30) g*(t )r] , is the number of relativistic degrees of freedom 
and rj(t ) ~ 10 4 — 10 5 is the ratio today of the energy density in matter to that in 
radiation. From eq. fl6.4Tj we can then see that accelerated expansion of the scale factor 
with a > is required to solve the horizon problem. As we have shown throughout 
the paper , such accelerated expansion can easily occur during the Hagedorn regime. 

We must here assume that, within an intitial horizon patch of size a p , there 
has been enough smoothing that it is sensible to talk about a single value of the 
temperature and entropy density. This same assumption must be made in every 
inflation model. 

6.2 The issue of non-adiabaticity 

Adiabatic Hagedorn inflating systems can in principle already solve the horizon prob- 
lem (as opposed to standard inflation in which nonadiabaticity is required). These 
systems could start with a large initial entropy due to the proximity of the initial 
temperature to the Hagedorn temperature. Subsequently, the entropy remains con- 
stant as the temperature drops and the scale factor a grows by a large amount. 
However in practice it is difficult to solve the horizon problem assuming adiabaticity: 
it is not possible to introduce a priori the entire entropy of the observable universe 
(i.e. S ~ 10 88 ) onto a primordial brane of order the string scale without either de- 
stroying the brane or having an unreasonably small string coupling (i.e. g s ~ 10 -88 !). 

Consider for example the L[— 1] system with a ~ t 4 / 3 . Then from the first entry 
in Table I, we see that the entropy on a world-volume a 3 is 

S p ~ (3 H a 3 p p ~ (3 H al J 1 (6.7) 

We will assume that the entropy on the 3-brane today is S oc (a to dayTtoday) 3 ■ 

We obtain the temperature/time relation during the Hagedorn phase from p ~ 
a~ 3 ~ i~ 4 by equating this expression for p with that from the first entry in Table I. 
We find tp = V^yffip - {3 H . Eq.Q can therefore be satisfied if 

f >^(V ±{ P P - M f i = ^". (6.8) 

'Jp J- today today 
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Hence even an adiabatic system here can in principle solve the horizon problem: 
provided that the primordial density is large enough (or the primordial temperature 
sufficiently close to Th) the horizon problem is solved even if So = S p . Indeed, we 
can see from eqn.( |6.7| ) that the initial entropy density can be as large as we like 
provided that we are prepared to accept a (3 that is arbitrarily close to (3h- 

The serious practical difficulty arises however when we consider the stability 
bound in eq. ( |4.14j ) since it clearly implies that non-perturbative effects will be im- 



portant unless the string coupling is fantastically weak. Indeed, generally if one 
assumes adiabaticity in attempting to solve the horizon problem, one requires the 
entire entropy of the observable universe to be present in string excitations at t = t p . 
Hence, if the volumes are initially of order the string scale, then the energy density 
must be enormous, p(t p ) ~ 10 88 (e.g. in eqn.|6~8l). After 60 e-folds of inflation p(t) 
falls below p c , the critical density which is needed to be in the Hagedorn phase, and 
the Yang- Mills phase takes over. However, for such a large initial value of the energy 



density, brane stability at t — t p (i.e. the constraint in eq.( 4.14 )) requires an initial 
value of g s < 1(T 88 . 

Because of this it is worth looking more critically at the supposed adiabaticity. 
This is a crucial assumption that almost certainly requires modification in a de Sitter 
background since the latter possesses a horizon with its own associated entropy. 
Hence, even if adiabaticity were the correct criterion (which it is not), it would only 
be meaningful for the coupled string/de Sitter system since there is a backreaction 
of the horizon on the string gas. The upshot for strings in a de Sitter background is 
that they are unstable to fluctuations and that this instability can sustain a period 
of de Sitter inflation |2TJ. This phenomenon is well known for strings once 



they are in the de Sitter-like phase 7 . However the missing ingredient that the present 
study adds is an explanation for how the universe enters a de Sitter like phase in 
the first place. Our mechanism gets the universe into the locally de Sitter phase, 
whereupon the mechanism of refs. |19|, ^TJ, keeps it there without having to assume 



adiabaticity. 

To achieve a sustained period of inflation therefore, we could just invoke the 



findings of refs. [19], ^(J, |2T|. However for the remainder of this section we briefly 
elaborate on this property of strings, using purely heuristic arguments, in order to 
indicate a possible direction for future study. 

Let us return to the density of states and consider instead a truly static metric 
in which the relevant Killing field is the time-translation, 

K^l. (6.9) 



7 We qualify 'de Sitter' because the exponential solutions we have do not possess the full 0(1,4) 
de Sitter symmetry. 
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De Sitter space is an example of a static metric; in an appropriate coordinate system 
the metric is 

ds 2 = -dt 2 (l - H 2 r 2 ) + TTT^dr 2 + rW (6.10) 

1 — H z r A 

and the Killing field has components (1, 0, 0, 0). We can therefore discuss thermody- 
namics by compactifying on an imaginary time coordinate with it = it + (3. 

A suitable modification to discuss the present case is to add an additional y 
coordinate with a metric g 55 = b(y) 2 which is time independent in this coordinate 
system and also independent of r; 

ds 2 = -(1 - H 2 r 2 )dt 2 + Ir^dr 2 + r 2 dtt 2 + biyfdy 2 . (6.11) 

1 — H 2 r 2 

Under this assumption each 4 dimensional leaf of the foliated space is a de Sitter space 
and the geometry unambiguously fixes the Hawking temperature of the Horizon, 
2n/H = (3. 

It is well established (see for example ref.|]22| and references therein) that the 
combined (geometric plus thermal) system obeys a generalized second law which is 
that the total entropy, 

S S ma tter ^horizon S ma tf er -\- ^ ^ -^horizon > (6-12) 

obeys 5S > 0, where A horizon is the area of the horizon. The constant-time surface 
S extends to the horizon r w 1/ H = f3 /2tc so that 

V\\ = Vol(S du )r d " ~ $\A honzon = VoliS^y^ 1 ~ ^n- 1 (6.13) 

where VollSd^) is the volume of the unit sphere in d\\ dimensions. 

Thus, assuming equilibrium at the temperature the total entropy of an L[-l] 
string gas in a de Sitter background is of the form 

s - v^-h) 2 + w ' (6 - 14) 

This function has a minimum at 

T crit = T H (l-0 ((£ 2 /^) 1/3 )) . (6.15) 

Below T crit fluctuations tend to drive the universe towards low temperature, i.e. out 
of the de Sitter phase 8 . The horizon has a negative specific heat which becomes 
relatively larger as the temperature drops so that equilibrium can be maintained 
(the condition being \CyJ\ > CV + ). Above T crit in the presence of a thermal gas 

8 We recognize that, in making this heuristic argument, we arc on rather thin ice since we are 
simultaneously assuming a string coupling which is small enough for our approximations to be 
correct, but large enough to maintain equilibrium between the different degrees of freedom. 
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of strings the flow is in the opposite direction towards an asymptotically limiting 
Hubble constant, H asymp = 2tiTh- Thus the generalized second law indicates that 
fluctuations in the geometry tend to drive the universe even further into the de 
Sitter phase provided that the Hubble constant is greater than the critical value. 
The entropy that the horizon loses when the Hubble constant is increased is more 
than offset by the huge increase in string entropy. Moreover, since the specific heat of 
the string gas becomes larger when the temperature increases, inevitably equilibrium 
is lost and energy flow to the strings becomes apparently limitless. 

There is a slight difference between this picture and the results of ref. which 
we should comment on. In ref.[20|, the density never goes above a critical density, 
p' c , and indeed asymptotes to it from below as one goes back in time. In our case, 
on the other hand, we must have p > p c ~ 1 for the calculation to be valid (i.e. to 
be in the Hagedorn regime). Indeed in the present paper the energy density diverges 
as we approach the Hagedorn temperature. The difference is because ref.|2l| was 
concerned with stretched strings and introduced a cut-off in the momentum integral 
in order to find the fractal dimension of the string. This 'coarse graining' put an 
artificial upper limit on the amount of string that can be packed into the volume. In 
fact they found the maximum fractal dimension (i.e. 2) only at the critical density p' c . 
A fractal dimension of 2 is typical of the random walk behaviour which exists once 
Hagedorn behaviour sets in. So the coarse graining in ref. [ 20fl effectively removed 
the Hagedorn behaviour and consequently ref. |2rJ| found that p = p' c gave (3 = (3h- 
Conversely, the present paper begins in the regime p > p c where the calculations of 
ref. PHI end. 



7. Conclusion and discussion 

In this paper we have studied the possible cosmological implications of the Hagedorn 
regime of open strings on D-branes in the weak coupling limit. Our main result in 
sections 2-5 is that, due to the non-extensive dependence of the free energy on the 
volumes, a gas of open strings can exhibit negative pressure leading naturally to a 
period of power law or even exponential inflation - Hagedorn inflation. We also find 
that the open string gas can dominate the cosmological evolution at weak coupling 
even though the D-brane tension becomes large in this limit. 

Hagedorn inflation also has a natural exit since any significant cooling can cause 
a change in the thermodynamics if winding modes become quenched or if the density 
drops below the critical density, p c ~ 1, needed for the entropy of the Hagedorn 
phase to be dominant. Such a cooling can be caused by a sudden adiabatic increase 
in the transverse radius or by the inflation itself. We find this 'easy-exit' feature 
of open-string Hagedorn inflation to be of its most appealing features. In addition, 
we found that a small but negative cosmological constant, can cause the universe to 
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enter a stable but oscillating phase. The effect of the Hagedorn phase is to soften 
the singular behaviour associated with the bounce. 

The most striking aspect of our discussion is probably the existence of negative 
pressure. One might therefore ask how general a feature this is expected to be. By 
T-dualizing we argued that we can put the negative pressure down to the fact that 
in any particular direction, the gravitational degrees of freedom have both Kaluza- 
Klein modes and winding modes whereas the open strings (which are dominant in 
the entropy) have only one or the other. Hence we expect negative pressure to be 
possible whenever there are large space-filling modes that dominate the entropy. 

In section 6 we speculated on how the inflation might be sustained through the 
well known phenomenon of string instability in a de Sitter background. 

Hagedorn inflation may be thought of as a first example in the search for alter- 
natives to the cosmological constant within the framework of string/brane systems. 
One possible direction for further study in this area is connected with the fact that 
we have throughout been taking the string coupling to be weak enough so that the 
brane tension does not play an important role in the cosmology. It is therefore in- 
teresting to ask if new cosmological effects might arise from large scale fluctuations 



in the branes themselves [y, pq] . The arguments of section 6 indicate that if it does 
then the brane driven inflation may be qualitatively different to the string driven 
inflation discussed here. This is because open strings are one dimensional objects 
with S ~ E + constVE whereas fluctuating p-branes have an entropy |24 



S~E^. (7.1) 

On calculating (3 = dS/dE we see that p ^ 1 branes do not have the divergent 
behaviour which is the defining feature of strings. The thermodynamics of these 
objects has been the subject of much study |24|] and there is probably more to be 
learnt here. 

Finally, an interesting connected issue which we did not discuss is related to the 
effect of brane melting discussed in ref. || , in which the non-perturbative aspect of 
D-brane thermal production must be taken into account. At the present time it is 
hard to make any quantitative estimate of these effects on the cosmological solutions 
we have been discussing here, but we hope to be able to address this question in 
future work. 
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2 




Figure 1: Examples of D-branes in thermal contact, showing (schematically) the direction 
of energy flow in order to achieve equilibrium, a: A non-limiting (open string excitations of 
D-brane) plus marginal limiting (closed string) system. Energy is lost into the transverse 
volume until the energy density on the non-limiting system falls below the critical Hagedorn 
density, b: A limiting (open string excitations of D-brane) plus marginal limiting (closed 
string) system. Energy flow is into open string winding modes as T — > Tjj. c: Limiting 
Dp brane plus limiting Dq brane plus marginal limiting (closed string) with p > q. Energy 
flows into open string windings on the Dp brane as T — > Th- 
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